feedback.t = 10;
feedback.rho = 0.98;
feedback.h = 300;
feedback.alpha = zeros(20,1);
feedback.beta  = zeros(20,1); 

settings.D_set   = 0.98;
settings.Fact_D  = 0.0;
settings.R2_set  = 40;
settings.Fact_r2 = 0.0;
settings.Ht_set  = 300;
settings.Fact_H  = 0.0;
settings.varR2_set = 200;
settings.Fact_varR2 = 1.0;
settings.P       = 5;
settings.m       = 20;
settings.dt      = 5;
%% Case #1
w1 = [0.2,0.2,0.2,0.2,0.2];
f1 = objfun103(w1,feedback,settings)

%% Case #2
% w2 = [0.1,0.15,0.2,0.25,0.3];
% f2 = objfun103(w2,feedback,settings);